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ABSTRACT 

We examine the dynamics and X-ray spectrum of the young Type la supernova remnant 0509—67.5 
in the context of the recent results obtained from the optical spectroscopy of its light echo. Our 
goal is to estimate the kinetic energy of the supernova explosion using Chandra and XMM-Newton 
observations of the supernova remnant, thus placing the birth event of 0509—67.5 in the sequence 
of dim to bright Type la supernovae. We base our analysis on a standard grid of one-dimensional 
delayed detonation explosion models, together with hydrodynamic and X-ray spectral calculations of 
the supernova remnant evolution. From the remnant dynamics and the properties of the O, Si, S, 
and Fe emission in its X-ray spectrum we conclude that 0509—67.5 was originated ~ 400 years ago 
by a bright, highly energetic Type la explosion similar to SN 1991T. Our best model has a kinetic 
energy of 1.4 x 10 51 erg and synthesizes 0.97 M© of 56 Ni. These results are in excellent agreement with 
the age estimate and spectroscopy from the light echo. We have thus established the first connection 
between a Type la supernova and its supernova remnant based on a detailed quantitative analysis of 
both objects. 

Subject headings: ISM: individual (SNR 0509—67.5) — nuclear reactions, nucleosynthesis, abundances 
— supernovae: general — hydrodynamics — supernova remnants — X-rays:ISM 



1. INTRODUCTION 

Astronomical observations can probe the material 
ejected by supernova (SN) explosions during two tran- 
sient phases with very different time scales. The initial 
optical transient (the SN itself) lasts for several months, 
and the ejecta structure can be studied through the emis- 
sion and absorption lines produced as the photosphere 
recedes into the exploded star. After the SN fades away, 
the ejected material starts to interact with its surround- 
ings, and the supernova remnant (SNR) phase begins. 
In this phase, the ejecta structure is revealed by the X- 
ray emission of the material heated by the reverse shock 
on timescales of hundreds or thousands of years. Both 
approaches are valid in principle, but up to now the dis- 
parity of the timescales involved has made it impossible 
to verify their mutual consistency by applying them to 
the same object. 

The serendipitous discovery of light echoes associated 
with the young SNRs 0509-67.5, 0519-69.0, and N103B 
in the Large Magellanic Cloud (LMC) by iRest et ail 
(2005) (henceforth R05) has opened new possibilities for 
establishing connections between SNRs and their parent 
SNe. In particular, spectroscopy of the light echoes has 
the potential to confirm the type of the SN explosion 
in a straightforward way, avoiding the difficulties inher- 
ent to typing SNRs from their X-ray spectra (for an in- 
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depth discussion of these diff iculties, see iRakowsk i et al.1 
120061 ). In a recent paper, IRest et alj 1)20071 ) (hence- 
forth R07) examined the light echo associated with SNR 
0509—67.5, and established that this object originated 
in a Type la explosion, in agreement with t he longstand- 
ing c l aims based on X-ray o bservations (|Hughes et al.l 
[19951 1 Warren fc Hughes! l200i henceforth WH041 . Fur- 
thermore, the quality of the light echo was such that it 
allowed for a meaningful comparison with the spectra of 
a variety of Type la SNe. Based on these comparisons, 
R07 concluded that the parent event of SNR 0509-67.5 
belonged to the group of overluminous, hi ghly energetic 
Type la SNe whose prototype is SN 1991T (jBenetti et all 
l2005h . 

These results offer an excellent opportunity to revisit 
the X-ray spectrum of SNR 0509-67.5. Recent de- 
velopments in modeling the thermal emission in young 
Type la SNRs using hydrodynamic calculations and 
no nequilibrium ionization simu lations (HD+NEI mod- 
els, iBadenes et al.|[2~003l . l2005al ) have led to a better un- 
derstanding of the relationship between the structure of 
SN ejecta and the X-ray spectra of SN Rs. By apply- 
ing HD+NEI models to the Tycho SNR, IBadenes et all 
( 2006) were able to estimate the kinetic energy and nu- 
cleosynthetic yields of the explosion. We can now do 
the same for 0509—67.5, comparing the results of the 
HD+NEI models with the light echo spectroscopy. In 
particular, we want to examine to what extent is it pos- 
sible to determine the brightness of a Type la SN by 
studying the X-ray emission from its SNR hundreds of 
years after the explosion. 

This paper is organized as follows. In § [21 we begin 
with a brief overview of delayed detonation models. In 
§[3[ we use these explosion models, together with the for- 
ward shock radius and velocity to constrain the dynam- 
ics of SNR 0509—67.5. In § [U we review the Chandra 
and XMM-Newton observations of 0509—67.5, and we 
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Fig. 1. — Nucleosynthetic yields of Fe (after radioactive decays, 
solid line, squares), Si (dotted line, triangles), and O (dashed line, 
diamonds) in the DDT models as a function of E^. Each model 
is labeled in the plot, indicating the corresponding value of ptr in 
units of 10 7 gcm -3 . 

calculate the values for the most relevant line centroids 
and line flux ratios. In § [5] we describe the HD+NEI 
models, and we evaluate them against both the dynamic 
and spectral constraints posed by the observations of 
0509—67.5. Finally, we discuss our results and present 
our conclusions in § [6l 

2. OVERVIEW OF DELAYED DETONATION MODELS 

Since they were first introduced by iKhokhlovl (|1991h , 
delayed detonations (DDTs) have become the most suc- 
cessful models for Type la SNe. In this kind of explo- 
sion, the burning front starts propagating as a subsonic 
flame (deflagration) in the central regions of a C+O white 
dwarf (WD) close to the Chandrasekhar mass, and then 
a transition to a supersonic regime (detonation) is arti- 
ficially induced, usually at a prescribed density p tr - The 
resulting ejecta structure is characterized by an approx- 
imately exponential density profile, with a composition 
dominated by Fe-peak nuclei in the center, surrounded 
by a shell rich in intermediate mass elements (IMEs: Si, 
S, Ar, Ca, etc.), and a thinner outer region dominated 
by O. Minor traces of C might remain at large radii, but 
most of the WD is burnt in the explosion. Spherically 
symmetric DDT models are able to explain the funda- 



mental properties of Type la SNe , including light curve 
shapes dHoflich fe Khokhfovill996l), optical and IR spec- 
tral evolution (IWheeler et al.l 119981: iBaron et alJ 120061 : 



iMarion et all l2006t iGerardv et al.ll2007ft. a nd spectropo- 
larimetric observations ( Wang et al. 20Q7|). 

The deflagration-to-detonation transition density p tr 
is the most important parameter in DDT explosions. 
Models with higher values of pt r are more energetic, 
and synthesize more Fc-peak elements, less IMEs, and 
less O than models with low values of pt r - Since Type 
la light curves are powered by the radioactive decay of 
56 Ni, more energetic m odels lead to more luminous SNe 
(jStritzinger et al.ll2006h . The luminosity of Typ e la SNe 
is als o tightly correlated to the light curve width ([Phillips 
|1993|) . although the physical processes connecting this 
parameter and the mass of 56 Ni synthesized in the ex- 
plosion are not as simp le as was initia lly th ought (see 
iKasen fe Wooslevl 120071: IWooslev etail [2007). In any 
case, the relationship between Type la SN luminosity 



(or light curve width) and the chemical st ructure of SN 
ejecta has now been firmly established by iMazzali et al.l 
(2007), and is consistent with the predictions of one- 
dimensional DDT models. 

The fundamental properties of the grid of DDT mod- 
els that we will use in the present work are plotted in 
Figure [Tl which can be com pared to the 'Zorro' dia- 
gram in lMazzali et all (120071). Models DDTa DDTc , 



and DDTe are taken from lBadenes et al.l (|2003l . l2005af l: 
mo del DDTg has been calculated with t he same code 
(see lBravo et al.lll996l ; IBadenes et al.ll2003f) . and extends 
the grid to lower values of Ek- The total masses of 
56 Ni synthesized in the models are 0.97 M Q (DDTa), 
0.74M Q (DDTc), 0.51 M (DDTe), 0.29 M Q (DDTg). 
To put these values in context, the estimat ed 56 Ni masses 
in the 20 objects considered 'normal' by IMazzali et al.l 
(2007) range between 0.94 ± 0.05 M (for SN 1991T) and 
0.24±0.05M o (for SN 1991M). A theoretical upper limit 
to the mass of radioactive Ni that can be obtained from 
the thermonuclear burning of a Chandrasekhar C+O 
WD is set by the p rompt detonation model DET in 
iBadenes et ail (|2003h . which yields 1.16 M of 56 Ni with 
E k = 1.6 x 10 51 erg. 

3. SNR DYNAMICS 

The dynamics of 0509—67.5 are constrained by two 
pieces of observational evidence. The first is the angu- 
lar radius measured by Chandra (15.1", see Table 3 in 
IBadenes et al.l 120071 ). which translates into a linear ra- 
dius of 1.1 x 10 19 cm (3.7 pc) at the known distance 
to the LMC (50 kpc. I Alvesll2004l ) . This measurement 
can be considered very accurate, because the errors on 
both angular radius and distance are extremely small (at 
the few percent level). The second observable is the for- 
ward shock (FS) velocity inferred by modeling the broad 
component of the Ly/3 e mission line, 3600 — 7100 kms -1 
(|Ghavamian et al.l f2007). Being a less direct measure- 
ment, the FS velocity is more subject to uncertainty 
than the FS radius, so this value should be considered 
wi th some measure of caution (see discussion in § 5.1 
of IBadenes et al.ll2007| ). To reproduce the dynamics of 
0509—67.5, we have used a one-dimensional hydrody- 
namic code to simulate the interaction between the ejecta 
density profiles from the DDT explosion models in our 
grid and a uniform ambient medium (AM; for more de- 
tails on the mo dels and a justificatio n of the uniform AM 
hypothesis, see Badene s et al.l [2007h . A successful SNR 
model must be able to match both FS radius and veloc- 
ity for a reasonable combination of age t and AM density 
RAM- 

The comparison between hydrodynamic models and 
observations, however, is not completely straightforward, 
because our simulations do not include the effect of cos- 
mic ray (CR) acceleration. This physical process can 
strongly modify the dynamics of the SNR, slowing down 
the FS and reducing the gap between FS and contact dis- 
continuity (CD) by an amount that dep ends on the ac- 
celeration efficiency (jEllison et al.ll2004h . Indeed, WH04 
found evidence for a significant nonthermal component 
in the X-ray spectrum of 0509—67.5 and argued that CR 
acceleration was taking place at the FS. Although the 
presence of nonthermal emission does not necessarily im- 
ply that the dynamics are CR-modified, we consider this 
to be a likely possibility. It is important to note that CR 
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Fig. 2. — Constraints on the dynamics of 0509—67.5 from the 
observed FS radius and velocity. The horizontally striped areas in- 
dicate the regions of the parameter space where the FS radius lies 
between Rfs and 1.1 X Rqd in models DDTa (blue) and DDTg 
(red). The vertically striped areas indicate the regions of the pa- 
rameter space where the FS velocity lies between upg and uqd 
in each model. The areas where both constraints overlap are high- 
lighted with solid colors. The age estimate from the light echoes is 
represented with the vertical solid and dotted lines. 

acceleration at the FS does not have an impact on the 
dynamics of the shocked ejecta, and in particular, it does 
not affect the position of the CD (at least to first order, 
see Fig. 2 in lEllison et all 120071) . This means that the 
CR-modified SNR radius must lie between the FS radius 
Rfs and the CD radius Rod in hydrodynamic models 
without CR acceleration. For one-dimensional models, 
it is necessary to allow for the fact that Rayleigh- Taylor 
instabilit ies at the CD increase th e projected radius by 
~ 10% (|Wang fe Chevafierl l200lh . Likewise, the CR- 
modified FS velocity must lie between the FS velocity 
ufs and the CD velocity ucd- We have summarized 
these constraints in Figure [2] for the most and least en- 
ergetic models in our DDT grid, DDTa and DDTg. The 
DDT models require that the SNR age be between 250 
and 610 yr to reproduce the FS dynamics, in very good 
agreement with the completely independent 400 ± 120 yr 
estimate from the light echoes (R05). The allowed val- 
ues of pam he between 5 x 10~ 26 and 2 x 10~ 24 gcm~ 3 , 
with more energetic models demanding higher AM den- 
sities. The FS radius provides the strongest dynamical 
constraints, imposing a tight correlation between pam 
and age that we shall revisit in § 15.31 

4. X-RAY OBSERVATIONS 

4.1. Chandra and XMM-Newton Data Sets 

SNR 0509-67.5 has been observed by the CCD cam- 
eras on both Chandra and XMM-Newton. The most de- 
tailed study of its X-ray emission to date was presented 
by WH04, who analyzed a Chandra ACIS-S data set 
taken on 2000 May 12-13 (ObsID 776, PI J.Hughes). The 
net exposure time was 47.9 ks, and the image was taken 
on the back-illuminated S3 chip. The authors found 
that the X-ray spectrum is dominated by line emission 
from the SN ejecta, with a small continuum contribution, 
probably nonthermal emission from the FS. The centroid 
of the Si Ka line blend indicates an unusually low (below 
He-like) ionization state for the plasma. The abundances 
obtained from plane-parallel shock fits favor a Type la 
origin, in agreement with the previous qualitative analy- 



sis of lHughes et al.l (|1995f) . These abundances were com- 
pared to_Hie_vjeUs of two type la SN explosion models 
from llwamoto et al.l ()1999f ). a 'fast deflagration' and a 
DDT, with the fitted values showing some preference for 
the DDT model, albeit with gross discrepancies in the Fe 
abundance that were noted by the authors. This might 
be due to the simplicity of the plane shock models used 
by WH04, but the poor statistics of the Fe Ka blend in 
the Chandra ACIS-S data set also make it difficult to 
constrain the properties of shocked Fe in 0509—67.5. 

For the present work, we have recalibrated the Chandra 
ACIS-S data using the latest versions of CIAO (3.4) and 
CALDB (3.4.0), which yielded 48.5 ks of useful exposure. 
To complement the Chandra data set, we have also re- 
duced and analyzed the XMM-Newton observation taken 
on 2000 July 4 (ObsID 0111130201, PI M.Watson). The 
instrument modes in this observation were small win- 
dow for EPIC-MOS1 (frame time: 300 ms), large win- 
dow for EPIC-MOS2 (900 ms) and EPIC-pn (48 ms). 
The medium filter was used. We reduced the data us- 
ing the latest SAS version (7.1.0) and calibration files. 
We applied the flare reject ion procedure described in 
iCassam-Chenai et al.l ((2004s) , which left a total exposure 
of 32.6 ks for all cameras. To create spectra, we selected 
single and double events (pattern < 4) for the EPIC-pn 
camera, which greatly improved the statistics in the Fe 
Ka line. The spectra were then rebinned to achieve a 
signal-to-noise ratio > 3o\ 

A preliminary inspection of the calibrated and re- 
duced spectra reveals that the centroids of the brightest 
lines in the EPIC-MOS cameras appear shifted towards 
high energies compared to both EPIC-pn and Chandra 
ACIS-S. The Si Ka centroid values obtained with the 
different CCD instruments, with standard 90% confi- 
dence intervals, are 1.850 ± 0.002 keV (EPIC-MOS1), 
1.847±o:oo2 keV (EPIC-MOS2), 1.834±0.002 keV (EPIC- 
pn), and 1.833±g;ggl keV (ACIS-S). The overlap between 
the ACIS-S and EPIC-pn values is important, because 
the energy scale of the Chandra ACIS-S data set was veri- 
fied using on-board calibration sources (see § 2 in WH04) . 
Furthermore, EPIC-pn has a substantially larger effec- 
tive area at the Fe Ka line than EPIC-MOS. In view of 
this, we have chosen to maximize the consistency of our 
data sets by not including the EPIC-MOS spectra in our 
analysis. The spatially integrated Chandra ACIS-S and 
XMM-Newton EPIC-pn spectra are plotted in Fig. [3] 

4.2. Spectral Fits to the Line Emission 

We have fitted the Chandra ACIS-S and XMM-Newton 
EPIC-pn data between 1.6 and 8.0 keV with a spectral 
model consisting of a power law continuum and 16 Gaus- 
sian lines, including K-shell transitions from principal 
quantum level n — 2 (Ka blends) for Si, S, Ar, Ca, and 
Fe; K-shell transitions from n — 3, 4, 5 (K/3, K7, and K(5) 
for Si, S, and Ar, and the Lya lines for the H-like ions of 
Si and S. Only the centroids and fluxes of the Ka blends 
and the Si K/3 line have been fitted. The centroids of the 
other K/3 lines and all the K7, K<5, and Lya lines have 
been kept fixed at the nominal energies for He-like ions. 
The fluxes of the K7 and KS lines have been tied to the 
respective K/3 lines taking the ratios at T e = 5 x 10 7 K 
from the ATOMDB data base ((Smith et al.ll2001h . The 
resulting fits are shown on Figure [4] and the fitted pa- 
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3. — Spatially integrated spectrum from Chandra ACIS-S 
(blue) and XMM- Newton EPIC-pn (red). The expected centroids 
of the Ka blend from the He-like ions of O, Ne, Mg, Si, S, Ar, and 
Ca are indicated for reference, as is the Fc Ka blend and the region 
dominated by Fe L-shell lines. 

rameters are listed in Table [TJ 

Outside the energy range of the fits, we have calculated 
a conservative upper limit to the O Ka flux by assum- 
ing that all the photons between 0.55 and 0.58 keV come 
from this blend, with interstellar absorption set at the 
maximum value found by WH04, N H — 0.076 x 10 22 cm 2 . 
We note that the presence of O was necessary in the 
plane shock fits performed by WH04, and that O lines are 
clearly seen in the grating observations (Hughes et al. in 
preparation), but for this element it is especially hard to 
disentangle the ejecta contribution from the shocked AM 
contribution (if any). Thus, we settle for the conservative 
upper limits listed in Table[T]to constrain our ejecta mod- 
els. We have also calculated the flux in the broad band 
from 0.8 keV to 1.2 keV, which is dominated by L-shell 
lines from Fe, with some contribution from K-shcll Nc 
lines, assuming a nominal value of Nh = 0.07 x 10 22 cm 2 . 
The principal line flux ratios normalized with respect to 
Si Ka are listed in Table H 

The values listed in Tables [T] and [2] merit a few com- 
ments. All the line centroids and line flux ratios, as 
well as the power law indices and normalizations, are 
consistent in both data sets within the statistical uncer- 
tainties. However, some puzzling differences between the 
Chandra ACIS-S and the XMM-Newton EPIC-pn spec- 
tra can be appreciated both in the Tables and in Fig- 
ures [3] and [U The fluxes of the three brightest lines are 
statistically inconsistent, with deviations of 11% for Si 
Ka, 24% for S Ka, and a very significant 34% for Si 
K/3. The upper limits for the O Ka emission and the 
fluxes in the broad Fe L band also show disagreements of 
9% and 18%. Some of these differences are larger than 
the cross-calibration uncertainty between Chandra and 
XMM-Newton, which is 15% in normalisation between 
0.8 and 2 keV (M.Stuhlinger et al., document XMM- 
SOC-CAL-TN-0052, Issue 3.0, January 2006). The dis- 
crepancies certainly deserve some further investigation, 
but they are of no consequence for the scientific objec- 
tives of this paper, for two reasons that will become clear 
along § [5] First, we only model parameters that are 
consistent in both data sets. And second, the tolerance 
ranges that we use to compare models and data are much 
larger than the differences between ACIS-S and EPIC-pn 



TABLE 1 

Spectral parameters in the SNR 0509—67.5 



Parameter 



Chandra XMM-Newton 
ACIS-S EPIC-pn 



Power Law Continuum 



a 

Norm. 



3.37 



+0.23 
0.21 



3.42^ 



(10- 6 photonscm- 2 s" 1 ) 4161*4 


348l«? 




Line Fluxes (10 6 phot cm 2 s 




OKa a 
Fe L b 
Si Ka 
Si K/3 
S Ka 
Ar Ka 
Ca Ka 
Fe Ka 


< 396 
1177 
123 ±4 
11.5 ± 1.8 
33.3 ±2.8 
2.94 ± 1.02 
0.992+ - 71 
4.601?:°° 


< 362 
958 
110 ±3 
7.64 ±0.89 

25 4+ 1 9 
3.56 ±0.88 
1.08 ±0.46 

q OQ + 0. 66 

°- oz -0.B5 


Line Centroids (keV) 


Si Ka 
Si K/3 
S Ka 
Ar Ka 
Ca Ka 
Fe Ka 


-, +0.003 

1.000_Q QQ2 

2 146+ 019 
z.i*o_ 014 

o 41 c+0.006 

3 064+ ' 027 

O.UD4_ 024 
„ _„_+0.072 
-0.060 

6.572_o.208 


1.834 ±0.002 
2.161 ±0.014 

9 499 +0.006 
z.i*zz_ Q05 

3.036 ±0.042 

01 1+0.032 
o.oll_ 

fi 440+0.030 

o.44u_ 02g 



Goodness-of-Fit 


X 2 /dof 


45.1/64 


135.0/130 



Note. — The limits given are the formal 90% confidence 
ranges (Ax 2 = 2.706). 

a Upper limit to the flux in the O Ka blend (0.55-0.58 keV) 
assuming N H = 0.076 X 10 22 cm~ 2b Flux in the 0.8-1.2 keV 
band assuming Nh — 0.07 x 10 22 cm -2 , including unresolved 
contributions from Ne Ka and other lines. 



TABLE 2 
Line flux ratios in the SNR 0509- 



67.5 





Chandra 


XMM-Newton 


Line Ratio 


ACIS-S 


EPIC-pn 



O Ka/Si Ka 
Fe L/Si Ka 
Si K/3/Si Ka 
S Ka/Si Ka 
Ar Ka/Si Ka 
Ca Ka/Si Ka 
Fe Ka/Si Ka 



< 3.2 
9.6 

0.093 ±0.018 

0.27 ± 0.03 
0.024 ±0.009 
0.0081 ± 0.0060 
0.037 ±0.017 



< 3.3 
8.7 

0.069 ±0.010 

0.23 ± 0.02 
0.032 ± 0.009 
0.0098 ± 0.0044 
0.030 ±0.007 



Note. — For simplicity, symmetric confidence ranges 
have been calculated taking the largest deviations from 
the fitted values in each case. 

that we have discussed here, and could easily accomodate 
the EPIC-MOS data sets as well. 

5. SPECTRAL MODELING 

5.1. Method, Parameters, and Strategy 

To compare our grid of DDT models to the observa- 
tions of SNR 0509—67.5, we have generated synthetic X- 
ray spectra for the shocked ejecta emission in our hydro- 
dynamic simulations followi ng the methods presented in 
iBadenes et al.l (120031 1200 5a). with updated atomic data 
( Badenes et al.l2006l) , and including radiative and ioniza- 
tion losses with an isobaric approximation as described 
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Fig. 4. — Spectral fits to the line emission in 0509-67.5 between 1.6 and 8.0 keV for Chandra ACIS-S (left) and XMM-Newton EPIC-pn 
(right). The expected centroids of the Ka blends from the He-like ions of Si, S, Ar, and Ca, and the position of the Fe Ka blend are 
indicated for reference. 
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Fig. 5. — Fundamental line flux ratios (left panels) and line centroids (right panels) as a function of pam in the DDT models at t = 400 
yr. For each model, the different values of f) are plotted with solid lines and diamonds (/3 m i n ), dotted lines and triangles (/3 = 0.01), and 
dashed lines and squares (/3 = 0.1). The line flux ratios are normalized to the observed values (from Tabled, so that 1.0 is a perfect 
match. Tolerance ranges for line ratios and centroids are indicated by the striped regions around the observed values (vertical solid lines). 
The horizontal dotted lines mark the values of pam that satisfy the dynamic constraints at t=400 yr, as discussed in § [3] 



in iBadenes et al.l ((2007) . For a given Type la SN explo- 
sion model, the synthetic spectra are controlled by three 
variables only: AM density pam, SNR age t, and j3. The 
parameter (3 represents the amount of collisionlcss elec- 
tron heating at the reverse shock, and is defined as the 
ratio of specific internal energies in electrons and ions at 
the shock transition. It can take values between (3 m in, 
which represents mass proportional heating, and 1, which 
represents full equilibration (for more details, see § 2.2 in 



IBadenes et al.ll20"05aD . We have performed simulations in 
a grid of seven pam values (5 x 10~ 26 , 10~ 25 , 2 x 10~ 25 , 
5 x 10~ 25 , 10~ 24 , 2 x 10~ 24 , and 5 x lO^gcm" 3 ) and 
three (3 values (/3 m m, 0.01, and 0.1). A similar, less ex- 
tended grid of synthetic X-ray spectra based on DDT 
models was presented and discussed in IBadenes et al.l 
(|2005bD . 

The ability of our synthetic spectra to reproduce SNR 
observations is limited by the quality of the atomic data 
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in the spectral code. In the particular case of 0509—67.5, 
several important issues arise due to the low ionization 
state of the plasma. At ionization states below He-like, 
the atomic data for the Ka blends of the major elements 
are reasonably complete, but deficiencies in K-shell tran- 
sitions from higher levels should be expected in all ele- 
ments. Atomic data in the Fe L complex are notably de- 
ficient, and altogether absent for ionization stages below 
p e +i6 (N e _fikB Fe). The strongest Fe L line in the Chan- 
dra HETG spectrum (at ^0.73 keV) can be associated 
with Ne-like Fe, and there are no signs of a significant 
contribution from lower ionization stages (Hughes et al., 
in preparation). Nevertheless, Fe L emission in the syn- 
thetic spectra should always be considered with caution 
as a matter of principle. 

Given the limitations of the synthetic spectra listed 
above and the issues with the X-ray data described in 
§ 14.21 we have chosen to focus our efforts on modeling 
only fundamental quantities for which we can trust both 
models and observations. These include the centroids 
of the three brightest line blends (Si Ka, S Ka, and Fe 
Ka), and the O Ka/Si Ka, Fe L/Si Ka and Fe Ka/Si 
Ka line flux ratios, noting again that special care must 
be taken in the case of the Fe L/Si Ka ratio. As we shall 
see, the combination of these parameters can constrain 
the kinetic energy of the SN explosion in the framework 
of the DDT models quite well. We will begin by using 
the DDT models at the nominal age of 400 yr in § 15.21 
and then discuss variations of age in § 15.31 and alternative 
models in § 15.41 In S I5.5l we will go beyond our selection of 
diagnostic parameters to evaluate the performance of the 
models across the entire spectral range of XMM-Newton. 

5.2. DDT Models at Nominal Age (t=400 yr) 

A comparison between the observed line flux ratios and 
centroids and the predictions of the synthetic spectra 
generated by the DDT models in our grid at t — 400 
yr is presented in Figure [5] In these plots, the observed 
values (vertical solid lines) are the averages of the Chan- 
dra ACIS-S and XMM-Newton EPIC-pn values, except 
in the case of the Fe Ka line, where only the XMM- 
Newton EPIC-pn data are taken into account. We note 
that the largest difference between ACIS-S and EPIC- 
pn in the relevant parameters corresponds to the the Fe 
L/Si Ka ratio, and is only 9%. Tolerance ranges around 
the observed values are highlighted with striped regions: 
a factor 2 above and below for line flux ratios (except 
for the O Ka/Si Ka ratio, where only the upper limit is 
constrained), 1% for the Si Ka and S Ka centroids, and 
0.5% for the Fe Ka centroid. These tolerance ranges are 
rather large, mostly because of systemati c uncertainties 
that are very hard to quantify (see § 5.2 in Badenes et al. 
120061 for a discussion), and will only be used to provide 
some guidance in the comparisons between models and 
data. For each DDT model, the values of pam that sat- 
isfy the dynamical constraints from § [3] at t = 400 yr are 
indicated by horizontal dotted lines. 

The combination of the O Ka/Si Ka and Fe Ka/Si 
Ka flux ratios clearly favors models with low O and high 
Fe emission. Only the DDT models with higher kinetic 
energies (and hence high Fe and low O yields), DDTa 
and DDTc, can reproduce both observables at the same 
time. When the dynamical constraints are brought into 
play, model DDTa (for pam = 10~ 24 g cm -3 and (3 be- 




FlG. 6. — Global view of the performance of the DDT models 
(mapped onto an axis) with varying pam at i = 400 yr. The 
striped regions represent the areas of the parameter space where 
the models can reproduce the observed O Ka/Si Ka (blue), Fe L/Si 
Ka (green), and Fe Ka/Si Ka (red) line flux ratios for an optimum 
value of p. The dynamical constraints on pam are overlaid in 
gray. The star indicates the position of model DDTa at pam = 
10 _24 gcm — 3 , where all the constraints overlap. 

tween 0.01 and 0.1) performs better than model DDTc, 
because of its lower O Ka flux. The Fe L/Si Ka ratio 
also favors model DDTa over DDTc, although we stress 
that comparisons based on this parameter must be of 
a qualitative sort. The high Fe L/Si Ka values found 
at low pam in models DDTc, DDTe and DDTg, for in- 
stance, are due mostly to Ne K emission, not to L-shell 
lines from shocked Fe. The line centroids plotted in the 
right hand panels of Figure [5] support the choice of model 
DDTa, but we note that the Fe Ka centroid is slightly 
underpredicted at the pam imposed by the dynamical 
constraints. 

The interplay between the constraints from the differ- 
ent line flux ratios with varying pam can be seen more 
clearly in Figure [6l where the sequence of DDT models 
has been mapped onto an £7. axis. The O Ka/Si Ka ra- 
tio can be reproduced by any model, provided that pam 
is high enough to ionize most O beyond the He- like stage 
at t = 400 yr. On the other hand, the Fe Ka/Si Ka ratio 
can only be reproduced by models with high Ek, where 
enough Fe has been shocked at t = 400 yr. These two flux 
ratios, together with the dynamical constraints clearly 
single out model DDTa with pam = 10~ 24 gcm~ 3 . The 
Fe L/Si Ka ratio confirms the choice, but we note again 
that the allowed region at low pam for this parameter is 
spurious. This plot highlights the remarkable agreement 
between the dynamical and X-ray spectral constraints on 
Pam for the energetic DDT models, specially consider- 
ing the large dynamic range of the diagnostic quantities 
shown in Figure This is not a trivial coincidence, and 
it is not observed in other explosion models, as we shall 
see in !5.4[ It indicates that the fundamental properties of 
the DDT models (exponential density profile, stratified 
ejecta) are well suited for SNR 0509—67.5. 

5.3. DDT Models at Different Ages 

The uncertainties in the age estimate from R05 (±120 
yr around the nominal age of 400 yr) stem from the un- 
known three-dimensional structure of the dust that is 
reflecting the light echo. The SNR dynamics are con- 
sistent with this range, but they impose a tight corre- 
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lation between t and pam, as shown in § [3j Now we 
shall see that the properties of the X-ray emission from 
0509—67.5 can substantially narrow down the age esti- 
mate, making large deviations from t — 400 yr seem un- 
likely. Higher ages require higher values of pam in order 
to reproduce the FS radius, and vice versa (see Figure[2]). 
These changes in pam lead to major differences in the 
ionization timescale r = J n e dt of the shocked ejecta, 
which rapidly drive several spectral parameters out of 
agreement with the observations. Two examples are pre- 
sented in Figure [3 where the O Ka/Si Ka flux ratio and 
the centroid of the Si Ka blend are compared with the 
observed values across the (pam , t) parameter space for 
model DDTa. At ages below 400 yr, most of the O is in 
the He-like ionization stage for the pam values imposed 
by the FS radius. This makes the O Ka emission too 
strong, with O Ka/Si Ka ratios in excess of 5 times the 
observed value. Above 400 yr, the O Ka/Si Ka ratio 
is in the allowed region, but then the He-like ion of Si 
becomes dominant in the plasma, pushing the centroid 
toward higher energies that are incompatible with the 
value measured by Chandra and XMM-Newton. These 
comparisons are less favorable for the other DDT mod- 
els, because their higher O content and lower Ek tend 
to increase the disagreement between the dynamical and 
spectral constraints. 

One last constraint on the age of 0509—67.5 comes from 
the lack of recorded historical evidence for a recent SN 
in the LMC. A normal to bright Type la SN at a dis- 
tance of 50 kpc would remain visible for several months, 
with an apparent visual magnitude around —1.0 at max- 
imum lig ht (the average extin ction towards the LMC is 
0.3 mag, llmara fc Blitz 2007). Such an event would be 
hard to miss by even the most inattentive observers. Eu- 
ropean exploration of the southern hemisphere was well 
under way 400 years ago, led by skilled navigators who 
relied heavily on astronomy and monitored the night sky 
whenever possible. Pieter Dirkszoon Keyser and Fred- 
erik de Houtman mapped the southern sky (including 
Mensa and Doradus, which encompass the LMC) from 
Java between 1595 and 1597, noting the positions of all 
the stars brighter than 5 th magnitude. In this island, the 
Dutch colony of Batavia (present day Jakarta, at 6° 16' S 
latitude) and the English colony of Bantam became per- 
manent around 1600. In South America, several Spanish 
settlements were thriving decades before that, including 
Lima (founded 1535, 12°12.6' S) and Buenos Aires (per- 
manent since 1580, 34°36' S), from which the LMC is 
circumpolar. Even allowing for some incompleteness in 
the historical records, a spectacular astronomical event 
like a bright LMC supernova should have left some kind 
of trace if it happened at any time later than the begin- 
ning of the 17 th century. 

Taken together, the light echoes from the SN, the dy- 
namic and spectral properties of the SNR, and the his- 
torical considerations, suggest an age very close to 400 yr 
for the birth event of 0509—67.5. The uncertainty around 
this value is hard to quantify, but deviations larger than 
a few decades seem unlikely, specially towards younger 
ages. 

5.4. Other Models 

We have seen that one-dimensional DDT models can 
reproduce the fundamental properties of the ejecta emis- 



sion in 0509—67.5 with remarkable accuracy. It is out- 
side the scope of this paper to perform an exhaus- 
tive exrjk2ration_of other Type la explosion paradigms 
like IBadenes et ail (|2006ft did for the Tycho SNR, but 
for the sake of completeness we discuss here some re- 
sults obtained using other models. In Figure [8] we 
present the line flux ratios at t = 400 yr predicted 
by a sub-Chandras ekhar explosion (model SCH from 
IBadenes et al.l l2003h an d a well-mixed 3D de flagration 
(model b30_3d_768 from iTravaglio et al.ll200l . In each 
case, we have overlaid the dynamical constraints on pam 
obtained with the procedure explained in § [3l The sub- 
Chandrasekhar model clearly shows that the overlap be- 
tween spectral and dynamical constraints found in the 
energetic DDT models is not trivial, and indeed does 
not happen for the ejecta structure of this edge-lit WD 
detonation. The well-mixed 3D deflagration model has 
more severe problems: the presence of Fe and O every- 
where in the ejecta leads to a systematic overprediction 
of the O Ka/Si Ka and Fe L/Si Ka flux ratios at all 
values of pam ■ 

5.5. Spatially Integrated Spectra 

In the previous Sections we have taken the approach 
of focusing on a few parameters that can be both re- 
liably extracted from the observations and confidently 
modeled with the existing spectral codes. This is a 
more meaningful and robust way of comparing synthetic 
spectra generated with HD+NEI models to data than 
the conventional method of \ 2 fi tting across the entire 
X-ray band (for dis cussions, see IBadenes et al.l [2005a; 
IBadenes et al.l 120061) . Nevertheless, it is always instruc- 
tive to study the performance of the spectral models 
from a more global point of view. In Figure [5J we plot 
the synthetic spectra of the DDT models at t = 400 
yr for the best value of pam within the dynamical con- 
straints, together with the XMM-Newton EPIC-pn data 
set. In each model, the value of (3 has been chosen to 
give the best approximation to the Fe Ka/Si Ka flux 
ratio: /3 = 0.02 for model DDTa and /3 = 0.1 for mod- 
els DDTc, DDTe, and DDTg. No spectral fitting of any 
kind has been done. Instead, the ejecta models have been 
normalized to match the Si Ka line and then a power 
law with the parameters determined in § 14.21 (a = 3.42, 
Norm = 3.5 x 10~ 4 photonscm~ 2 s _1 ) has been added 
for the continuum. Absorption has been set to the fidu- 
cial Njj — 0.07 x 10 22 cm 2 in all cases. By foregoing any 
spectral fitting, we can more easily compare the proper- 
ties of the ejecta models to the data and to each other. 

The performance of our model of choice, DDTa 
[pAM — 10~ 24 gcm -3 , (3 — 0.02) is very good at all en- 
ergies, albeit with some problems in localized areas. The 
flux in the Fe L complex is somewhat low, and the shape 
shows deviations from the observed spectrum. These dis- 
crepancies might be due mostly to the limitations in the 
atomic data discussed in ij l5.ll Other problems are prob- 
ably related to the explosion model, such as the deficit 
in the Mg Ka flux and the low centroid of the Ca Ka 
blend 5 . The degradation of the synthetic spectra along 
the sequence of DDT models with decreasing Ek (or p tr ) 

5 The same two issues were noted by Badcncs et al. (2006) in 
comparisons between the synthetic spectrum from model DDTc 
(t = 430 yr, p AM = 2 X 10- 24 gcm" 3 , /3 = 0.03) and the X-ray 
emission from the Tycho SNR. This suggests that the Mg content 
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Fig. 7. — Comparison between the observations and the values predicted by the synthetic spectra produced from model DDTa for the 
O Ka/Si Ka flux ratio (left panel) and the centroid of the Si Ka blend (right panel) as a function of pam and t. The black dashed plots 
are the dynamical constraints from the SNR radius from Figure [2] The age estimate from the light echoes is represented with the vertical 
solid and dotted lines. 

The light echo observations of R05 and R07 have 
turned 0509—67.5 into a unique object: the only X-ray 
bright SNR whose Type la origin is confirmed and whose 
explosion energy can be estimated based on spectroscopic 
information from its parent supernova. Motivated by 
these groundbreaking observations, we have re-examined 
the dynamics and X-ray spectrum of 0509—67.5, ap- 
pl ying the HD+NEI m odelin g techniques i ntrodu ced 
in IBadenes etUI ((2003f t and IBadenes et"all (|2005aj ) to 
Chandra and XMM- Newton observations. We have 
based our analysis on a grid of one-dimensional DDT 
explosions, taking advantage of the relationship between 
ejecta structure and in these models to place the birth 
event of 0509—67.5 in the sequence of dim to bright Type 
la SNe. Our conclusions are in excellent agreement with 
the light echo results of R07: SNR 0509-67.5 was orig- 
inated by a bright, highly energetic Type la event of 
the subtype often referred to as SN 1991T-like. We thus 
present, for the first time, a scenario where the dynamics 
and X-ray spectrum of a Type la SNR form a consistent 
picture with the spectroscopy of its parent supernova. 

Our preferred Type la explosion model, DDTa, has 
a kinetic energy of 1.40 x 10 51 erg, and synthesizes 
1.03 M of Fe, 0.09 M of Si, 0.07 M of S, 0.04 M 
of O, and O.14M of other elements (mostly Ar and 
Ca). The peak bolometric luminosity of this model 
is logLboi(ergs s _1 ) = 43.39, and the amount of 
56 Ni in the ejecta before nuclear decays is 0.97 M . 
These values compare very well with the estimates fo r 
SN 1991T: logL bol = 43.36 (iContardo et all [2000); 
M 56Ni = 0.87 M rr, dStritzinger et all 120061 UVOIR 
meth od), 0.94 M (|Stritzinger et alj|2006t iMazzali et all 
I2007L B-band lightcurve method). The ejecta density and 
chemical composition profiles in model DDTa (pam = 

-24 
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Fig. 8. — Fundamental line flux ratios as a function of pam 
form models SCH (top panels) and b30_3d_768 (bottom panels) at 
t = 400 yr. All labels and plots are as in Figure [5] 

is plain to see in Figure O specially in the behavior of the 
Fe and O emission. Normalization distances can be cal- 
culated for each model, yielding 36 kpc (DDTa), 51 kpc 
(DDTc), 23 kpc (DDTe), and 29 kpc (DDTg). Given 
the coarseness of our pam grid, values within ~ 50% of 
the distance to the LMC can be considered satisfactory. 
A more accurate study of the normalization distances 
cannot be justified without the inclusion of multidimen- 
sional effects in the HD+NEI simulations (see § 8.1 in 
IBadenes et aT1l2006l ). 

We conclude this Section with a brief mention of the 
continuum emission. Detailed modeling of this compo- 
nent of the X-ray spectrum is outside the scope of the 
present work, but our HD+NEI calculations can help 
clarify the issue of its origin raised by WH04. At high 
energies (> 2 keV), the flux of the thermal AM emission 
in our HD+NEI models (for pam = 10~ 24 gcm -3 and 
t = 400 yr at a distance of 50 kpc) is a factor 6 below 
the power law models listed in Table [TJ This is a strong 
indication that the origin of the continuum emission in 
0509—67.5 is indeed nonthermal. 

6. DISCUSSION AND CONCLUSIONS 



in our DDT models might be underestimated, and the distribution 
of Ca in the ejecta might be biased towards low densities. 



10 gem , /3 = 0.02) can reproduce very well both the 
dynamics (FS radius and velocity) and the fundamental 
properties of the X-ray spectrum (O, Si, S, and Fe emis- 
sion) of 0509—67.5 at an age of 400 yr. Several dynami- 
cal, spectral, and historical arguments indicate that the 
SNR age cannot be very different from this value, which 
also coincides with the completely independent estimate 
derived by R05 from the geometry of the light echo. 

These resul ts, together with our previous work on the 
Tycho SNR (jBadenes et al.l I2006D . constitute fir m evi- 
dence that the phenomenological one-dimensional DDT 
models, which have proved so successful in explaining the 
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Comparison between the XMM-Newton EPIC-pn data set and the DDT model spectra at t = 400 yr for the most favorable 
and pam (in units of 10 _24 g cm -3 ) within the dynamical constraints. Each ejecta model has been normalized to match the 
Si Ka blend, and then a power law model with the parameters listed in Table [T] has been added. Absorption has been fixed at 
X 10 22 cm 2 . The position of the principal emission lines (Ka and Lya for O, Ka for all other elements) has been indicated for 



Fig. 9. 
values of /3 
flux in the 
N H = 0.07 
clarity. 

light curves and spectra of Type la SNe, are equally ca- 
pable of reproducing the fundamental properties of the 
X-ray emission from young Type la SNRs. In particular, 
it is possible to use HD+NEI simulations based on these 
DDT models to estimate the kinetic energy (and ulti- 
mately the brightness) of a Type la SN explosion from 
the X-ray spectrum of its SNR. Hundreds of years af- 
ter the SN explosion, the memory of the cataclysmic 
event persists in the X-rays from its SNR, opening a 
window onto ages past. This enables us to explore the 
relationship between individual dim/bright Type la SNe 
and their immediate surroundings (presence or absence 
of stellar formation, local metallicities, etc.) with much 
grea ter detail than is available to extragalactic studies 
(see lPrieto et al.ll2007t ). In the context of SNR research, 
having a good estimate for Ek can help to build better 
models for the impact of CR acceleration on the SNR 
dynamics. HD+NEI simulations have clearly become a 
powerful and flexible tool to study the relationship be- 
tween SNRs and their parent SNe. For interested readers, 
the synthetic X-ray spectra from our models are available 
from the authors upon request. 

We conclude with a reminder that much is left to do 
in the study of Type la SNe and their SNRs in general, 
and of 0509—67.5 in particular. Ongoing direct mea- 
surements of the ejecta and FS expansion using grat- 
ing observations and proper motion studies should im- 
prove our knowledge of the dynamics of this SNR in 
the near future. Further work on the X-ray emission 
should take the asymmetries of the object into account, 



in particular the enhanced Fe knots found by WH04. 
The one-dimensional models presented here can only 
sample the average or bulk properties of the shocked 
eject a. Observational evidence for modera te asymme- 
tries J^se^^Fal. 2007; Gerardv et al.ll2007l) and clump- 
ing ((Leonard et a.l.ll2005t IWang et al.ll2007l ) in the ejecta 
of Type la SNe is widespread, and SNR studies offer 
a unique opportunity to study these effects. This might 
provide crucial constraints for the multi-dimensional sim- 
ulations of physical ly-motivated DDT explosions cur- 
rently under way (iBravo &: Garria-Send I2006I : iPlewal 
[2007t [Jordan et al.ll2007n ~ 
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Fig. 10. — 'The Persistence of Memory', by Salvador Dali, 1931. Oil on canvas, Museum of Modern Art, New York City. This painting 
by the genial Catalan artist represents the passage of time, and is very relevant to our work on SNR 0509—67.5 in a convoluted, surreal 
way. Dali referred to the vaguely human figure in the center as the 'paranoiac-critical Camembert', for whatever that is worth. 

APPENDIX 

SOME REFLECTIONS ON THE RELEVANCE OF SURREALISM TO ASTROPHYSICS 

The meaning of Salvador Dah's painting The Persistence of Memory (Figure I10[) has been discussed at length 
elsewhere (see f.i. MoMA Highlights, New York: The Museum of Modern Art). The most common interpretation is 
that the famous 'soft clocks' are a metaphor for the passage of time, and the ants represent the inevitable consequence 
thereof: corruption of the flesh. Other interpretations are definitely possible, and equally valid - this being, after 
all, surrealism. Dali once commented that the idea of the soft clocks came to him from the observation of a piece 
of Camembert cheese melting in the summer heat. Be that as it may, if the soft clocks do represent the passage of 
time, it is easy to draw an analogy between The Persistence of Memory and young SNRs like 0509—67.5. The X-ray 
observations of these objects open a window onto ages past, telling us what processes went on during the supernova 
explosions that originated them, and further back in time, how the supernova progenitor modified its surroundings. 
Just like a piece of melting Camembert cheese. 



